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The effective-action formalism is applied to a gas of bosons. The equations describing the con- 
densate and the excitations are obtained using the loop expansion for the effective action. For a 
homogeneous gas the Beliaev expansion in terms of the diluteness parameter is identified in terms 
of the loop expansion. The loop expansion and the limits of validity of the well-known Bogoliubov 
and Popov equations are examined analytically for a homogeneous dilute Bose gas and numerically 
for a gas trapped in a harmonic-oscillator potential. The expansion to one-loop order, and hence 
the Bogoliubov equation, is shown to be valid for the zero-temperature trapped gas as long as the 
characteristic length of the trapping potential exceeds the s-wave scattering length. 
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\ The dilute Bose gas has been subject to extensive study for more than half a century, originally in an attempt to 
(y-j . understand liquid Helium II, but also as an interesting many-body system in its own right. In 1947, Bogoliubov showed 
how to describe Bose-Einstein condensation as a state of broken symmetry, in which the expectation values of the field 
operators are non- vanishing due to the single-particle state of lowest energy being macroscopically occupied, i.e., the 
K>" ' annihilation and creation operators for the lowest-energy mode can be treated as c- numbers |^ . In modern terminology, 
' the expectation value of the field operator is the order parameter and describes the density of the condensed bosons. In 
Bogoliubov's treatment, the physical quantities were expanded in the diluteness parameter Vriof?, where tiq denotes 
Q.^ the density of bosons occupying the lowest single-particle energy state, and a is the s-wave scattering length, and 
Bogoliubov's theory is therefore only valid for homogeneous dilute Bose gases. The inhomogeneous Bose gas was 
^sj studied by Gross and Pitaevskii who independently derived a nonlinear equation determining the condensate 
. density. A field-theoretic diagrammatic treatment was applied by Beliaev to the zero-temperature homogeneous 
■ dilute Bose gas, showing how to go beyond Bogoliubov's approximation in a systematic expansion in the diluteness 
^ , parameter v^noo^ ^ ; and also showing how repeated scattering leads to a renormalization of the interaction between 
^ the bosons. This renormalization was in Beliaev's treatment a cumbersome issue, where diagrams expressed in terms 
of the propagator for the noninteracting particles are intermixed with diagrams where the propagator contains the 
interaction potential. Beliaev's diagrammatic scheme was extended to finite temperatures by Popov and Faddeev 
Q and was subsequently employed to extend the Bogoliubov theory to finite ternperatures by incorporating terms 

. containing the excited-state operators to lowest order in the interaction potential |Q, ^ . 
^ ' A surge of interest in the dilute Bose gas due to the experimental creation of gaseous Bose-Einstein condensates 
, occurred in the mid-nineties The atomic condensates in the experiments are confined in external potentials, which 
' poses new theoretical challenges; especially, the Beliaev expansion in the diluteness parameter \/riQa^ is questionable 
5h when the density is inhomogeneous. The renormalization of the potential was generalized to a trapped system by 
Proukakis et al. ||l^, [Tl|]. Leading-order corrections to the Gross-Pitaevskii equation for a trapped Bose gas were 
studied by Stenholm |12| . The finite temperature Beliaev-Popov theory was applied to a trapped gas by Fedichev 
and Shlyapnikov JlSf . 

In this paper we shall employ the two-particle irreducible effective-action approach, and show that it provides an 
efficient systematic scheme for dealing with both homogeneous Bose gases and trapped Bose gases. We show how the 
effective-action formalism can be used to derive the equations of motion for the dilute Bose gas, and more importantly 
that the loop expansion can be used to determine the limits of validity of approximations to the exact equations of 
motion in the trapped case. 

The paper is organized as follows. The model and the two-particle irreducible effective-action approach are intro- 



duced in Section fll The homogeneous dilute Bose gas is considered in Sec. Ill, and familiar equations of motion are 
rederived. In Sec. we demonstrate how the renormalization of the interaction potential due to repeated scattering 
is conveniently carried out in the effective-action formalism. In Sec. ^ we review the trapped Bose gas. The main 



results of the paper are presented in Sec. VI, where the equations of motion are solved numerically in order to assess 
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the limits of validity of approximations to the exact equations of motion. Finally, in Sec. VI] we summarize and 
conclude. 

II. EFFECTIVE ACTION FORMALISM FOR BOSONS 

We consider a system of spinless bosons described by the action 

S[tP, = J drdt^^r, t) [idt - i/o(r) + ^A ^(r, t) 

- \j drdr'dt'ilj\r,t)'iP\r',t)U{r~r')'4;{r',tyj{r,t) (1) 

where ip is the scalar field describing the bosons. Here fi denotes the chemical potential, Ho = p^/2m + ^(r) is the 
one-particle Hamiltonian consisting of the kinetic term and an external potential, and U{r) is the potential describing 
the interaction between the bosons. We have chosen units so that h = 1, but will restore Ti in final results. It will 
prove convenient to introduce a matrix notation whereby the field and its complex conjugate are combined into a 
two-component field 4> = {ip , tp^ ) = (i^i, 02)- 

The correlation functions of the Bose field are obtained from the generating functional 

Z[ri, K] = J V<j>exp (^tS[(l)] + + ^(f>^ K(P^ , (2) 

by differentiating with respect to the source = {rj,ri*) = {rji, 772). In Eq. (^, matrix notation is implied in order to 
suppress the integrations over space and time variables. A two-particle source term, K, has been added to the action 
in the generating functional in order to obtain equations involving the two-point Green's function in a two-particle 
irreducible fashion. 

The generator of the connected Green's functions is 

W[fj,K] = -tlnZ[r^,Ki (3) 

and the derivative 

5W 

gives the average field, 0, with respect to the action S'[0] -I- ri''(j) + (t)''K(j)/2, 

4>{r,t) = ^ = y l?00(r,i)exp (^zS[4,] + cP + '-cP^ Kc^^ = {cP{r,t)). (5) 

The average field <& is seen to specify the condensate density and is referred to as the condensate wave function. 
The derivative of W with respect to the two-particle source is 



+ ^G,,(r,i,r',t'), (6) 
where G is the full connected two-point matrix Green's function describing the bosons not in the condensate, 
r (rfr' ( mr,t)Si,Hr',t')) {S^{r,t)6i^{r' ,t')) \ 

where (SV'(r, t) is the deviation of the field from its mean value. Sip — tjj ~ $. Likewise, we shall write (p = (p + S(p for 
the two-component field. We note that in the path integral representation, averages over fields, such as in Eq. (j^), 
are automatically time-ordered. 

We introduce the effective action, F, the generator of the two-particle irreducible vertex functions, through the 
Legendre transform of the generator of connected Green's functions, W: 

F[<^,G] = W[fi,K]-7j^^-^^^K4>-'-TTGK. (8) 
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The effective action satisfies the equations ST/Scj) = —77 — K(f> and ST/SG = —iK/2. In a physical state where the 
external sources vanish, rj = = K, the variations of the effective action with respect to the field averages <f> and G 
vanish, yielding the equations of motion: 

^£.0, (9.) 

|§=0. (Sb) 

According to Cornwall et al. , the effective action can be written on the form 

r[0,G] = S[^] + ^TrlnGoG-i + '-Tr{G^' - E«)G - ^Trl + Tsi^, G], (10) 
where Go is the noninteracting matrix Green's function, 

Go-(r, t, r', ^ - ( f " + ^ - Ho + M ) - ^'^'^^ " ^^^^ 



and the matrix 



(50t(r,f),5(/)(r',t') 



+ Go-i(r,t,r',i'), (12) 



will turn out to be the self-energy to one- loop order (see Eq. (|2l|)). Using the action describing the bosons, Eq. (j^), 
we obtain for the components Il\^\r,t,r' ,t') = 6(t — t')I^[^\r,r'), where 



4\\r 


r') 


= ^(r' 


I]«(r 


r') 


= U{r 


I]«(r 


r') 


= U{r 


^22' 


r') 


= 6{r' 



) y dr"C/(r-r")|$(r",<)|2 + J7(r-r')$*(r,<)$(r',t). (13) 

The delta function in the time coordinates reflects the fact that the interaction is instantaneous. Finally, the quantity 
r2 in Eq. (^ is 

Ea = -Un(e'^-'["^''^"*l)^Pi, (14) 

where Sint[4>, S^i] denotes the part of the action S[(l) + 50] which is higher than second order in S(f) in an expansion 
around the average field. The quantity F2 is conveniently described in terms of the diagrams generated by the action 
•Sinti^, <^'/'], and consists of all the two-particle irreducible vacuum diagrams as indicated by the superscript "2PI", 
and the diagrams will therefore contain two or more loops. The subscript indicates that propagator lines represent 
the full Green's function G, i.e., the brackets with subscript G denote the average 



*^'"'['^'*'^1)g = (detiG)-i/2 / V{6(j))ei^'''^^ '54,^^S„A^M] _ (15) 



The diagrammatic expansion of F2 corresponding to the action for the bosons, Eq. (|l]), is illustrated in Figure |i| 
where the two- and three-loop vacuum diagrams are shown. Since matrix indices are suppressed the diagrams are 
to be understood as follows. Full lines represent Green's functions and in the cases where we display the different 
components explicitly, Gn will carry one arrow (G22 can according to Eq. (|^) be expressed in terms of Gn and thus 
needs no special symbol), G12 has two arrows pointing inward and G21 carries two arrows pointing outward. Dashed 
lines represent the condensate wave function and can also be decorated with arrows, directed out from the vertex to 
represent or directed towards the vertex representing The dots where four lines meet are interaction vertices, 
i.e., they represent the interaction potential U (which in other contexts will be represented by a wiggly line). When all 
possibilities for the indices are exhausted, subject to the condition that each vertex has two ingoing and two outgoing 
particle lines, we have represented all the terms of r2 to a given loop order. Finally, the expression corresponding to 
each vacuum diagram should be multiplied by the factor i""^, where s is the number of loops the diagram contains. 
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FIG. 1: Two-loop (upper row) and three-loop vacuum diagrams (lower row) contributing to the effective action. 



It is well-known that the expansion of the effective action in loop orders is an expansion in Planck's constant h . 
The first term S[(f>] on the right-hand side of Eq. ( |l0|) is referred to as the zero- loop term and the terms where the 
trace is written explicitly as one-loop terms, because they are proportional to and H^, respectively. We note that 
the presented effective action approach is capable of describing arbitrary states, including non-equilibrium situations 
where the external potential depends on time. Although we in the present paper shall limit ourselves to study a 
Bose gas at zero temperature the theory is straightforwardly generalized to finite temperatures. For example, in the 
Schwinger-Keldysh technique for treating general non-equilibrium states, the field is just attributed an additional 
index. For an account of treating arbitrary states we refer to Ref. ||l^; see also Ref. ?? for an application of the 
Schwinger-Keldysh technique to the dilute Bose gas. The equations of motion (^-9b) together with the expression 
for the effective action, Eq. (Q, form the basis for our subsequent calculations. 



III. HOMOGENEOUS BOSE GAS 



We shall now consider the case of a homogeneous Bose gas in equilibrium. The theory of such a system is well 
known, but the effective action formalism will prove to be a simple and efficient tool which permits one to derive the 
well-known results with particular ease. For the case of a homogeneous Bose gas in equilibrium, the general theory 
presented in the previous section simplifies considerably. The single-particle Hamiltonian Hq is then simply equal to 
the kinetic term, Ho{p) = p^/2m = Sp, and the condensate wave function $(r, t) is a time and coordinate-independent 
constant whose value is denoted by y/no, so that uq denotes the condensate density. The first term in the effective 
action, Eq. (^, is then 

Sm = i^lno - ^Uonl) J drdtl, (16) 

where Uq = J drU{r) is the zero- momentum component of the interaction potential. For a constant value of the 
condensate wave function, <&(r,t) = ^/no, Eq. (|l3|) yields 

^^'-[noUp no{Uo + Up))- ^''> 
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Varying, in accordance with Eq. (pa|), the effective action, Eq. (pij|), with respect to Uq yields the equation for the 
chemical potential 



fi = noUo + - 



[{Uo + Up)iGn{p) + G22{p)) + UpiGuip) + G^iip))] - (18) 



where the notation for the four-momentum, p — (j),lu), has been introduced. The first term on the right-hand side is 
the zero- loop result, which depends only on the condensate fraction of the bosons. The second term on the right-hand 
side is the one- loop term which takes the noncondensate into account. The term involving the anomalous Green's 



functions G12 and G21 will shortly be absorbed by the renormalization of the interaction potential (see Sec. IV). From 
the last term originate the higher loop terms which will be dealt with at the end of the section. 

The equation determining the Green's function is obtained by varying the effective action with respect to the matrix 



Green's function G{p), in accordance with Eq. (9b), yielding 

ST _ i 
dG~~2 

where 



0=-^ = -^f-G-i+Go-i+E«+I]'), (19) 



s:.-2z^. (20) 



Introducing the notation for the matrix self-energy E — E'^^ -|- E', Eq. ( |l9| ) is seen to be the Dyson equation 



G-^=Go^-Y.. (21) 

In the context of the dilute Bose gas, this equation is referred to as the Dyson-Beliaev equation. 

The Green's function in momentum space is obtained by simply inverting the 2x2 matrix Gq^{p) — E(p) resulting 
in the following components: 



Gii(p) 

Gl2(p) 
G2l(p) 



- £p - ^ + E22(p) 



Dp 



Dp ' 

-^2l{p) 



G22{P) = -- + ^r>-J^ + ^.dp) ^ (22) 



p 



all having the common denominator 

Dp = (cc; + ep-Ai + E22(p))(to-ep-f/i-Sn(p)) + 5]i2(p)S2i(p). (23) 

From the expression for the matrix Green's function, Eq. (|^), it follows that in the homogeneous case its components 
obey the relationships G22{p) = Gii(— p) and Gi2(— p) = Gi2(p) = G2i(p). The corresponding relations hold for the 
self-energy components. We note that the results found for fi and G to zero- and one-loop order coincide with those 
found in Ref. Q to zeroth and first order in the diluteness parameter -v/noo^. For example, according to Eq. (^7|) we 
obtain for the components of the matrix Green's function to one loop order 



uj^ -el- 2noC/pep 



G[\\P) .2 .2 



g1^2^(p)^ , 7%% J ^ (24) 
- - 2no(7pep 

which are the same expressions as the ones in Ref. [Q. As we shortly demonstrate, the loop expansion for the case of a 
homogeneous Bose gas is in fact equivalent to an expansion in the diluteness parameter. From Eq. (p^ ) we obtain for 

the single-particle excitation energies to one-loop order Ep — \J s^, + 2noC^p£p, which are the well-known Bogoliubov 

energies ffl. 
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Differentiating with respect to no the terms in r2 corresponding to the two-loop vacuum diagrams gives the two- 
loop contribution to the chemical potential. Functionally differentiating the same terms with respect to Gji gives 
the two-loop contributions to the self-energies Sy . The diagrams we thus obtain for the chemical potential and 
the self-energy E are topologically identical to those found by Beliaev however, the interpretation differs in that 
the propagator in the vacuum diagrams of Fig. ^ is the exact propagator, whereas in reference the propagator to 
one-loop order appears. 

In order to establish that the loop expansion for a homogeneous Bose gas is an expansion in the diluteness parameter 
Vnoo^, we examine the general structure of the vacuum diagrams comprised by Any diagram of a given loop order 
differs from any diagram in the preceding loop order by an extra four-momentum integration, the condensate density 
no to some power k, the interaction potential U to the power fc -I- 1, and k + 2 additional Green's functions in the 
integrand. We can estimate the contribution from these terms as follows. The Green's functions are approximated by 
the one-loop result Eq. (p^). The additional frequency integration over a product of k + 2 Green's functions yields k + 2 
factors of noJ7 (where U denotes the typical magnitude of the Fourier transform of the interaction potential), divided by 
2k + 3 factors of the Bogoliubov energy E. The range of the momentum integration provided by the Green's functions 
is (TOnoC/)i/2. The remaining three-momentum integration therefore gives a factor of order ng '^"''^^^rn'^/^J/ ^=+1/2^ 
and provided the Green's functions make the integral converge, the contribution from an additional loop is of the 
order {nom^U^y^'^ . This is the case except for the so-called ladder diagrams, in which case the convergence need to 
be provided by the momentum dependence of the potential. The ladder diagrams will be dealt with separately in the 
next section where we show that they through a renormalization of the interaction potential lead to the appearance of 
the t-matrix which in the dilute limit is proportional to the s-wave scattering length a and inversely proportional to 
the boson mass. The renormalization of the interaction potential will therefore not change the estimates performed 
above, but only change the expansion parameter. Anticipating this change we conclude that the expansion parameter 
governing the loop expansion is for a homogeneous Bose gas indeed identical to Bogoliubov's diluteness parameter 
\/noa^. 



IV. RENORMALIZATION OF THE INTERACTION 



Instead of having the interaction potential appear explicitly in diagrams, one should work in the skeleton diagram- 
matic representation where diagrams are summed so that the four-point vertex appears instead of the interaction 
potential, thus accounting for the repeated scattering of the bosons. In the dilute limit, where the inter-particle 
distance is large compared to the s-wave scattering length, the so-called ladder diagrams give the largest contribution 
to the four-point vertex function [ pT| . The ladder diagrams are depicted in Fig. ^ On computing the corresponding 
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FIG. 2: Summing all diagrams of the "ladder" type results in the t-matrix, which to lowest order in the diluteness parameter 
is a momentum independent constant g, diagrammatically represented by a circle. 



integrals, it is found that an extra "rung" in a ladder contributes with a factor proportional not to V nom^U^ as 
was the case for the type of extra loops considered in the previous section, but to komU, where fco is the upper 
cut-off momentum (or inverse spatial range) of the potential, as first noted by Beliaev The quantity kornU is not 
necessarily small for the atomic gases under consideration here. Hence, all vacuum diagrams which differ only in the 
number of ladder rungs that they contain are of the same order in the diluteness parameter, and we have to perform 
a summation over this infinite class of diagrams. The ladder resummation results in an effective potential T(p,p', q), 
which is called the t-matrix and is a function of the two ingoing momenta and the four-momentum transfer. Due 
to the instantaneous nature of the interactions, the t-matrix does not depend on the frequency components of the 
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ingoing four-momenta, but for notational convenience we display the dependence as T{p^p',q). To lowest order in 
the diluteness parameter, the t-matrix is independent of four- momenta and proportional to the constant scattering 
amplitude, T(0, 0,0) — Airh^a/m — g, where a is the s-wave scattering length This is illustrated in Fig. ||, 

where we have chosen an open circle to represent g. Iterating the equation for the ladder diagrams we obtain the 
well-known t-matrix equation: 

T{p,p\q) = U^ + l I dS'U^,Gii{p + q')Gn{p-q')T{p + q\p^q',q^q'). (25) 

If the theory is generalized to finite temperatures, the t-matrix takes into account the effects of thermal population 
of the excited states. 

We shall now show how the ladder resummation alters the diagrammatic representation of the chemical potential 
and the self-energy (see also Ref. [0). In Fig. ^ is displayed some of the terms up to two-loop order contributing to 
the chemical potential /x. The first two terms in Eq. ( p^ is represented by diagrams a-d, and the two-loop diagrams 
e-f originate from The diagrams labeled e and f are formally one loop order higher than c and d, but they differ 
only by containing one additional ladder rung. Hence, the diagrams c, d, e, and f, and all the diagrams that can be 





FIG. 3: Diagrams up to two-loop order contributing to the chemical potential. Only the two-loop diagrams relevant to the 
resummation of the ladder diagrams are displayed. The two-loop diagrams not displayed are topologically identical to those 
shown, but differ in the direction of arrows or the presence of anomalous instead of normal propagators. 



constructed from these by adding ladder rungs, are of the same order in the diluteness parameter y/noo^ as showed 
at the end of the previous section. They are therefore resummed, and as discussed this leads to the replacement of 
the interaction potential U by the t-matrix. 

We note that no ladder counterparts to the diagrams a and b in Fig. |^ appear explicitly in the expansion of the 
chemical potential, since such diagrams are two-particle reducible and are by construction excluded from the effective 
action F. However, diagram b contains implicitly the ladder contribution to diagram a. In order to establish this 
we first simplify the notation by denoting by Np the numerator of the exact normal Green's function Gii(p), which 
according to Eq. ( p^ ) is Np = uj + Sp — ^ + Sii(— p). We then have Dp — NpN-p — Ei2(p)S2i(p) — D^p, and the 
contribution from diagram b can be rewritten on the forms 



/ d^pUpGi2{p) = / d^pUp ^h(pL 

J ^ J ^ P7VpiV_p-Si2(p)S2l(p) 



^MNpN^p _ Sl2(p)S2l(p)Sl2(p) 



= j d^Up (Ei2(p)Gii(p)Gii(-p) - ^2i{p)Gi2{pf) 

= J d^pUp {noUpGn{p)Gni-p) + [^Mp) - noUp]Gn{p)Gni~p) - ^2i{p)Gi2ipf) ■ (26) 








V 
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FIG. 4: Diagrammatic representation of the last two rewritings in Eq. ( |2q ) which lead to the conclusion that the diagram b 
of Fig. 1^ implicitly contains the ladder contribution to diagram a. The anomalous self-energy E12 is represented by an oval 
with two ingoing lines, E21 is represented by an oval with two outgoing lines, and the sum of the second and higher-order 
contributions to E12 is represented by an oval with the label "2" . 



In Fig. ^ the last two rewritings are depicted diagrammatically. We see immediately that the first term on the right- 
hand side corresponds to the first ladder contribution to diagram a, and since to one-loop order, Ei2(p) = noUp^ the 
other terms in Eq. (26) are of two- and higher- loop order. The self-energy in the second term on the right-hand side 
can be expanded to second loop order, and by iteration this yields all the ladder terms, and the remainder can be 
kept track of analogously to the way in which it is done in Eq. (p6|). The resulting ladder resummed diagrammatic 
expression for the chemical potential, displayed in Fig. ||, is seen to be equal to that found in Ref. p7[ |. 



-\i= O +2 




+ 



FIG. 5: The chemical potential to one-loop order after the ladder summation has been performed and the resulting t-matrix 
been replaced by its expression in the dilute limit, the constant g. 



In the same manner, the self-energies are resummed. For En a straightforward ladder resummation of all terms is 
possible, while for E12 the same procedure as the one used for diagrams a and b in Fig. || for the chemical potential 
has to be performed. In Fig. ^, we show the resulting ladder resummed diagrams for the self-energies En and E12 to 
two-loop order in the dilute limit where T{p,p' , q) w g. 

In Ref. a diagrammatic expansion in the potential was performed, which yields to first order the diagram 
E(f ) in Fig. " 



but not the other two-loop diagrams. This theory, where the normal self-energy is taken to be 

Ell = E^"'' + E^""*, the anomalous self-energy to E12 — E$^2°\ ^-^d the diagrams displayed in Fig. ^ are kept in the 
expansion of the chemical potential, is referred to as the Popov approximation. Although we showed at the end of Sec. 
Ill that all the two-loop diagrams of Fig. are of the same order of magnitude in the dilutcness parameter \/noa^ at 
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la 2a 2b 2c 2d 



FIG. 6: Normal, En, and anomalous, E12, self-energies to two-loop order after the ladder summation has been performed. 

zero temperature, the Popov approximation applied at finite temperatures is justified, when the temperature is high 
enough, kT 3> gnQ. Below, we shall investigate the limits of validity at zero temperature of the Popov approximation 
in the trapped case. 

In this and the preceding section we have shown how the well-known expressions for the self-energies and chemical 
potential for a homogeneous dilute Bose gas are conveniently obtained using the effective action formalism, where 
they simply correspond to working to a particular order in the loop expansion of the effective action. We have 
established that an expansion in the diluteness parameter is equivalent to an expansion of the effective action in the 
number of loops. Furthermore, the method provided a way of performing a systematic expansion, and the results are 
easily generalized to finite temperatures. We now turn to show that the effective action approach provides a way of 
performing a systematic expansion even in the case of an inhomogeneous Bose gas. 

V. INHOMOGENEOUS BOSE GAS 

We now review the experimentally relevant case of a Bose gas trapped in an external static potential, in order to 
set the stage for the numerical calculations in Section 0. In this case the Bose gas will be spatially inhomogeneous. 
The effective action formalism is equally capable of dealing with the inhomogeneous gas, in which case all quantities 
are conveniently expressed in configuration space as presented in section || . Varying, in accordance with Eq. (|9^), 
the effective action F, Eq. (p^), with respect to $*(r,f), we obtain the equation of motion for the condensate wave 
function: 

Q oFo 

{lUdt -Ho + M)$(r, t) - g |$(r, t)f <i>(r, t) + 2igGn{r, t, r, i)$(r, t) - . . (27) 

o<l>*(r, t) 

To zero-loop order, where only the first term on the right-hand side appears, the equation is the time-dependent Gross- 
Pitaevskii equation [|| ^ . We have already, as elaborated in the previous section, performed the ladder summation by 
which the potential is renormalized and the t-matrix appears and substituted its lowest order approximation in the 
diluteness parameter, the constant g. Since the t-matrix in the momentum variables is a constant in the dilute limit, 
it becomes in configuration space a product of three delta functions, r(ri, r2, ra, = gS{ri — r4)(5(r2 — r4)(5(r3 — r4). 
The quantity r2 is defined as the effective action obtained from F2 by summing the ladder terms whereby U is replaced 
by the t-matrix, and its diagrammatic expansion is topologically of two-loop and higher order. 

The Dyson-Beliaev equation, Eq. (^, and the equation determining the condensate wave function, Eq. (|2^), 
form a set of coupled integro-differential self-consistency equations for the condensate wave function and the Green's 
function, with the self-energy specified in terms of the Green's function through Eq. (^0|). The Green's function can 
be conveniently expanded in the amplitudes of the elementary excitations. We write the Dyson-Beliaev equation, Eq. 
on the form 



dr"dt" [ihdt<T3S{r - r")S{t - t") + (T3L(r, t, r", t")] G(r", i", r', t') = ;il,5(r - r')5{t - t'), (28) 
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where we have introduced the matrix operator L{r, t, r', t') = a3Ho5{r — r')5{t — t') + (T3E(r, t, r', t') and (73 is a PauH 
matrix. Up to one-loop order, the matrix S is diagonal in the time and space coordinates and we can factor out the 
delta functions and write L(r, t, r'. t') = S(t — t')S{r — r')L(r), where 

^^''>- [-g^*{v)^ -Ho + fi-2gmr)\^ )■ ^^■^> 

The eigenvalue equation for L are the Bogoliubov equations |l|. The Bogoliubov operator L is not Hermitian, but 
the operator a^L is, which renders the eigenvectors of L the following properties |2^, For each eigenvector 
(Pj{r) — {uj{r),Vj{r)) of L with eigenvalue Ej, there exists an eigenvector (pj{r) = {v* {r),u*{r)) with eigenvalue —Ej. 
Assuming the Bose gas is in its ground state, the normalization of the positive-eigenvalue eigenvectors can be chosen 
to be {ipj, ifk) = Sjk where we have introduced the inner product 



{^3^^k) ^ J drip](r)a3ipk{r) = J dr{u*{r)uk{r) - v*{r)vk{r)). (30) 

It follows that the inner product of the negative-eigenvalue eigenvectors (p are 

{^j,^k) = J <ir^](r)cr3^fe(r) = J dr{vj{r)vlir) - Uj{r)ulir)) = -Sjk (31) 

and the eigenvectors tp and tp are mutually orthogonal, {(pj,(pk) ~ 0. By virtue of the Gross-Pitaevskii equation, the 
vector (po{r) = (<i>(r), — $*(r)) is an eigenvector of the Bogoliubov operator L with zero eigenvalue and zero norm. 
In order to obtain a completeness relation, we must also introduce the vector (/3a(r) = ($a(r), — ^^(r)) satisfying the 
relation Ltpa — ctfo, where a is a constant determined by normalization, (cpo, fa) = 1- The resolution of the identity 
then becomes 

J2' {pAMir') - ^,(r)^](r')) a, + (^a(r)^S(r') + ^o(r)^i(r')) a, = lS{r r'), (32) 

j 

where the prime on the summation sign indicates that the zero-eigenvalue mode ipo is excluded from the sum. Using 
the resolution of the identity, Eq. (|3^), allows us to invert Eq. (p8|) to obtain the well-known spectral representation 
of the Green's function 

G(r,r',.) . nj:' (z^^.W^Kr') ~ ■ (33) 

It follows from the spectral representation of the Green's function, that the eigenvalues Ej are the elementary exci- 
tation energies of the condensed gas (here constructed exphcitly to one-loop order). Using Eq. (^), we can at zero 
temperature express the noncondensate density or the depletion of the condensate, rinc = n — uq, in terms of the 
Bogoliubov amplitudes 

n„c(r) = z j ^Gn(r,r,..) = Y^\v,{v)\\ (34) 



VI. LOOP EXPANSION FOR A TRAPPED BOSE GAS 



We now turn to determine the validity criteria for the equations obtained to various orders in the loop expansion 
for the ground state of a Bose gas trapped in an isotropic harmonic potential V{r) — ^muj'^r'^. To this end, we shall 
numerically compute the self-energy diagrams to different orders in the loop expansion. 

Working consistently to one-loop order, we need only employ Eq. (|2^) to zero loop order, providing the condensate 
wave function which upon insertion into Eq. (^9|) yields the Bogoliubov operator L to one-loop order, from which 
the Green's function to one-loop order is obtained from Eq. (p3|). The resulting Green's function is then used 
to calculate the various self-energy terms numerically. In order to do so, we make the equations dimensionless 
with the transformations r = aosc^, = ^/NoJa^^, Uj = Uo^J'^Uj, Ej = hutEj, and g = {hu;ta^g^/NQ)g, where 
flosc = y'n/mujt is the oscillator length of the harmonic trap, and No is the number of bosons in the condensate. 

To zero- loop order, the time- independent Gross-Pitaevskii equation on dimensionless form reads 



11 



We solve Eq. ( pSD numerically with the steepest-descent method, which has proven to be sufficient for solving the 
present equation p3[ |. The result thus obtained for $ is inserted into the one- loop expression for the Bogoliubov 
operator L, Eq. (|29|), in order to calculate the Bogoliubov amplitudes uj and vj and the eigenenergies Ej. Since 
the condensate wave function for the ground-state, $, is real and rotationally symmetric, the amplitudes Uj,Vj in 
the Bogoliubov equations can be labeled by the two angular momentum quantum numbers I and m, and a radial 
quantum number n, and we write Unim{f , 9 , (j)) = Uni{f)Yim{9 , (f) , Vnim{r,0,4>) = Vniir)YimiO, (j)). The resulting 
Bogoliubov equations are linear and one-dimensional: 

II d\ , 11(1 + 1) ,1.2 ~ , ^~^2rA~ , ~A2.~X- 



o - fi~2'' + o ~2 +0^-1^ + 25$ (^^) u^iir) + g1>\r)vni(r) = E„iu„i(r), 

1 1 52 ^ , 11(1 + 1) , 1^.2 . , o^A2.-x \ , ^^2,^.~. 



2 f gf2-^ + 2 f2 + 2^ ^^"' + 2.9$ (^)j = -EM(r). (36) 

Note that the only parameter in the problem is the dimensionless coupling parameter 5 = AnNoa/aosc- Solving the 
Bogoliubov equations reduces to diagonalizing the band diagonal 2M x 2M matrix L, where M is the size of the 
numerical grid. The value of M in our computations was varied between 180 and 360, higher values for stronger 
coupling, and the grid constant has been chosen to O.OSflosc giving a maximum system size of 18aosc- We have used 
MATLAB to perform the diagonalization. 

In the following we shall estimate the orders of magnitude and the parameter dependence of the different two- 
and three- loop self-energy diagrams, and to this end we shall use the one-loop results for the amplitudes u, v and 
the eigenenergies E obtained numerically. When working to two- and three-loop order, one must also consider the 
corresponding corrections to the approximate t-matrix g. These contributions have been studied in Ref. [pj]| , and their 
inclusion will not lead to any qualitative changes of our results. In fact, even at finite temperature, the dependence of 
the t-matrix on the coupling is weak as long as the temperature is not close to the critical temperature for Bose-Einstein 
condensation Tc 0, |l^. 

Let us first compare the one-loop and two-loop contributions to the normal self-energy. The only one-loop term is 

E[f\r, r', uj) = 25|$(r)|25(r - r') - 2gno(r)d(r - r'). (37) 

We first compare S^^"' with the two-loop term which is proportional to a delta function, i.e., the diagram 2a in Fig. 
^. We shall shortly compare this diagram to the other two- loop diagrams. For diagram 2a we have 

Slf)(r,r',c.) = 2*.g<5(r-r') J ^G(r,r,u;') = 2gn„e(r)5(r - r'). (38) 

The ratio of the two-loop to one-loop self-energy contributions at the point r is thus equal to the fractional depletion 
of the condensate at that point. In Figure ^is shown the numerically computed dimensionless fractional depletion at 
the origin, Ti,ic(0)/rio(0)j where we have introduced the dimensionless notation 

no(f) = |$(f)P, 

n,e(f)=E'l*^(f)l'- (39) 

j 

We have chosen to evaluate the densities at the origin, r = 0, in order to avoid a prohibitively large summation over 
the I 7^ eigenvectors. As apparent from Figure 0, the log-log curve has a slight bend initially, but becomes almost 
straight for coupling strengths g > 1000. A logarithmic fit to the straight portion of the curve gives the relation 

^nc(O) ^ ^ nn,n~1.2 



"o(0) 



0.00195^-^ (40) 



For weaker coupling, the curvature is seen to be stronger. When we reintroduce dimensions, the power-law relationship, 
Eq. (|4^), is multiplied by the inverse of the number of bosons in the condensate, N^^, because the actual and 
dimensionless self-energies are related according to 

where s denotes the loop order in question. The ratio between different loop orders of the self-energy is thus not 
determined solely by the dimensionless coupling parameter g = AnNoa/aosc, but by A'o and a/aosc separately. We 
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10° 10^ 10^ 10^ 

g 

FIG. 7; Fractional depletion of the condensate at the trap center as a function of the dimensionless coupling strength g — 
'inNoa/ aosc- Asterisks denote our numerical results, circles denote the local-density approximation with the numerically 
computed condensate density inserted, and the line is the local-density approximation using the Thomas-Fermi approximation 
for the condensate density. 



thus obtain for the fractional depletion in the strong-coupling limit, g ^ 1000, 

1.2 



nnc(O) _ 1 "-nc(O) 
no(0) ~ TVo no(0) 



0.041A^" 



a 

Q-osc 



(42) 



It is of interest to compare our numerical results with approximate analytical results such as those obtained using 
the local density approximation (LDA). The LDA amounts to substituting a coordinate-dependent condensate density 
in the expressions valid for the homogeneous gas. The homogeneous-gas result for the fractional depletion is [0 



no 30? 

In the strong coupling limit we can use the Thomas- Fermi approximation for the condensate density 

1 



»^o(r) 



/ISTVoay/' 

\ ^OSC / 



1 



ISiVoa 



2/5 p 



(43) 



(44) 



which is obtained by neglecting the kinetic term in the Gross-Pitaevskii equation (24[ . For the fractional depletion at 
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the origin there results in the local density approximation 

,(0) (15iVo)^/' f ^ V^' 



no(0) 3tt^V2 \a. 



(45) 



as first obtained in Rcf. |25j . The LDA is a valid approximation when the gas locally resembles that of a homogeneous 
system, i.e., when the condensate wave function changes little on the scale of the coherence length ^, which according 
to the Gross-Pitaevskii equation is f = (87r?io(0)a)^^/^. For a trapped cloud of bosons in the ground state, its radius 
R determines the rate of change of the density profile. Since i? is a factor g^/^ larger than ^ j2^, we expect the 
agreement between the LDA and the exact results to be best in the strong-coupling regime. The fractional depletion 
of the condensate at the trap center as a function of the dimcnsionless coupling strength g = A-k N^a / aosc is shown in 
Fig. ^ In Fig. ^ are displayed both the local-density result Eq. (^3|) with the numerically computed condensate density 
inserted, and the Thomas- Fermi approximation (|4g), showing that the LDA indeed is valid when the coupling is strong. 
Furthermore, inspection of Eq. (|45) reveals that the LDA coefficient and exponent agree with the numerically found 
result of Eq. (|4|), which is valid for strong coupling. However, when g < 10, the LDA prediction for the depletion 
deviates significantly from the numerically computed depletion. Inserting the numerically obtained condensate density 
into the LDA instead of the Thomas- Fermi approximation is seen not to substantially improve the result. 

The relation for the fractional depletion, Eq. (|4^), is in agreement with the results of Ref. ||l2|, where the leading- 
order corrections to the Gross-Pitaevskii equation were considered in the one-particle irreducible effective-action 
formalism, employing physical assumptions about the relevant length scales in the problem. These leading-order 
corrections were found to have the same power-law dependence on Nq and a/aosc- A direct comparison of the 
prcfactors cannot be made, because the objective of Ref. ||l^ was to estimate the higher-loop correction terms to the 
Gross-Pitaevskii equation and not to the self-energy. 

The two-loop term Y^y^^ can at zero temperature be ignored as long as n^c ^ ^Oi which is true in a wide, 
experimentally relevant parameter regime. The result for the fractional depletion Eq. ( |4^ ) depends very weakly on 
A^O) so as long as A'o does not exceed 10^, which is usually fulfilled in experiments, we can restate the criterion for 
the validity of Eq. (^2|) into the condition a <C aosc- In experiments on atomic rubidium and sodium condensates, 
this condition is fulfilled, except in the instances where Feshbach resonances are used to enhance the scattering length 
MM. 



In section III we showed that for a homogeneous gas all two-loop diagrams are equally important in the sense that 



they are all of the same order in the diluteness parameter ^/n^a? . The situation in a trapped system is not so clear, 
since the density is not constant. We shall therefore compare the five normal self-energy diagrams S^^° in Fig. 
^, to see whether they display the same parameter dependence and whether any of the terms can be neglected. In 
particular, the Popov approximation contains the diagram E^^"'' but neglects all other second-order diagrams, and 
we will now determine its limits of validity at zero temperature. Since diagram 2a contains a delta function, we shall 
integrate over one of the spatial arguments of the self-energy terms and keep the other one fixed at the origin, r = 0. 
We denote by i?^-'^ the ratio between the integrated self-energy terms j and 2a, 

/drSfi")(0,r,c^ = 0) 

In Fig. H, we display the ratios R^^^ for the different integrated self-energy contributions corresponding to the diagrams 
where j represents 2b and 2c. The contributions from diagrams 2d and 2e are equal and turn out to be exactly equal 
within our numerical precision to that from diagram 2c. Furthermore, inspection of the diagrams in Fig. ^ reveals 
that when the condensate wave function is real, the anomalous contribution E^2°^ is equal to Ti^if^] '^'rP and Ex2^'' 
are equal to S^^'^-', and T^^i-f^ is equal to t!-^-^\ In the parameter regime displayed here, the contribution from diagram 
2a is larger than the others by approximately a factor of ten, and displays only a weak dependence on the coupling 
strength. In the weak-coupling limit, ^ 1, it is seen that the terms corresponding to diagrams 2b-e can be neglected 
as in the Popov approximation, with an error in the self-energy of a few per cent. When the coupling gets stronger, 
this correction becomes more important. A power-law fit to the ratio R^'^^^ in the regime where the log-log curve is 
straight yields the dependence 

^ 0.028g^-i^ (47) 

which is equal to 0.5 when 'g « 10^; for g greater than this value, the Popov approximation is seen not to be valid. 
If the ratio between the oscillator length and the scattering length is equal to one hundred, Cosc = 100a, the Popov 
approximation deviates markedly from the two-loop result when Nq exceeds 10^, which is often the case experimentally. 
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g 

FIG. 8: Ratio between different dimensionless two-loop self-energy terms as functions of the dimensionless coupling strength 
g = 4:TvNoa/aosc- Asterisks denote the ratio i?'^''' as defined in Eq. ( ^6| ) and circles denote the ratio R^^'^\ The terms i?'^**' and 
are equal and turn out to be similar in magnitude to R^^''\ and are not displayed. 



In order to investigate the importance of higher-order terms in the loop expansion, we proceed to study the three- 
loop self-energy diagrams. We have found the number of summations over Bogoliubov levels to be prohibitively 
large for most three-loop terms; however, we have been able to compute the two diagrams S^^"^ and E^l^^ displayed 
in Fig. 1^, for the case where one of the spatial arguments is placed at the origin thereby avoiding a summation 
over I ^ components. We compare the diagrams E^^"' and 212°' the two-loop diagrams. As we have seen, 
diagrams S^^''-' , S^^^'^'' and E^i'^^ in Fig. ^ are similar in magnitude and dependence on g, as are the anomalous two- 
loop diagrams S^2" "^6 have therefore chosen to evaluate only diagrams S^^'''' and Sj'^"''. The results for the ratios 
'Sf^\o,r,uj = 0)/T,\f''^ (0, r, u = 0) and E^2''-' (0, r, = 0)/'s''^^\o,r,u! = 0), evaluated for different choices of r, are 
shown in Fig. nO. A linear fit to the log-log plot gives for the normal terms the coefficient 0.016 and the exponent 
0.76 when r = O.Soosc and the coefficient 0.0029 and the exponent 0.78 when r — Qosc, and for the anomalous terms 
with the choice r — Oosc the coefficient is 0.0015 and the exponent 0.82. Restoring dimensions according to Eq. ( |4l| ) 
we obtain 

11!^"'°°^^'"^°^ - 0-15A^o-"" ( —) • (48) 



Eir(0,aosc,c^ = 0) 



The ratio between three- and two-loop self-energy terms in the homogeneous case was in Sec. [II found to be pro- 
portional to ynoQ^- A straightforward application of the LDA, inserting the central density no(0) into this formula, 
yields the dependence cx NQ-'^{a/aosc)^''^ ■ This is not in accordance with the numerical result, Eq. (|4S 



15 




FIG. 9: Self-energy diagrams to three loop order which are evaluated numerically. 

although the self-energies were evaluated at spatial points close to the trap center. The discrepancy between the LDA 
and the numerical result is attributed to the fact that we fixed the spatial points in units of Oosc while varying the 
coupling g, although the physical situation at the point r — Oosc (and r — ^Qosc and r = |aosc respectively) varies 
when g is varied. It is possible that the agreement with the LDA had been better if the length scales had been fixed 
in units of the actual cloud radius (as given by the Thomas-Fermi approximation) rather than the oscillator length. 
However, the present calculation agrees fairly well with the LDA as long as the number of atoms in the condensate 
lies within reasonable bounds. Since iVo > 1 in the condensed state, Eq. (^) yields that E^^"'' ^ E^^^-* whenever 
the s-wave scattering length is much smaller than the trap length. We conclude that only when this condition is not 
fulfilled is it necessary to study diagrams of three-loop order and beyond. 

VII. CONCLUSION 

We have applied the two-particle irreducible effective- action approach to a condensed Bose gas, and shown that it 
allows for a convenient and systematic derivation of the equations of motion both in the homogeneous and trapped 
case. We have chosen to work at zero temperature, but the formalism is with equal ease capable of dealing with 
systems at finite temperatures and general non-equilibrium states. Beliaev's diagrammatic expansion in the diluteness 
parameter and the t-matrix equations are expediently arrived at with the aid of the effective-action formalism. We 
have shown that the parameter characterizing the loop expansion for a homogeneous Bose gas is equal to the diluteness 
parameter, the ratio of the s-wave scattering length and the inter-particle spacing. For an isotropic, three-dimensional 
harmonic-oscillator trap at zero temperature, the small parameter governing the loop expansion has been found to 
depend on the ratio between the s-wave scattering length and the oscillator length of the trapping potential, and has 
a weak dependence on the number of particles. The expansion to one-loop order, and hence the Bogoliubov equation, 
is found to provide a valid description for the zero-temperature trapped gas when the oscillator length exceeds the 
s-wave scattering length. We have compared our numerical results with the local-density approximation, which is 
found to be valid when the number of particles is large compared to the ratio between the oscillator length and the 
s-wave scattering length. Furthermore, we have found that all the self-energy terms of two-loop order are not equally 
large for the case of a trapped system: in the limit when the number of particles is not large compared to the ratio 
between the oscillator length and the s-wave scattering length, the Popov approximation has been shown to be a valid 
approximation. 

Recent experiments have probed the regime of large scattering lengths in rubidium and sodium condensates ??. In 
the limit of very strong interactions, the loop expansion itself is questionable and other theoretical methods must be 
invoked p^ . However, in the intermediate regime of moderately dense clouds, the self-energy corrections considered 
in this paper should become important and are possible to study experimentally. 
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FIG. 10: Ratio of dimensionless three-loop to two-loop self-energy diagrams as a function of the dimensionless coupling strength 
g = AnNoa/ttoac- Asterisks denote the ratio of the normal self-energy terms Nq'E^^-^^ fE^^^ ^ evaluated at the point (0, Oosc, w = 
0), open circles denote the same ratio evaluated at (0, 0.5aosc, w = 0), and diamonds denote the same ratio evaluated at 
(0, l.SffloscW = 0). Crosses denote the ratio of anomalous self-energy terms A^'oE^j^V^i?' (0,aoBc, w = 0). 
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